On the intermediate superconducting state in a flat slab 
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The intermediate state (IS) in a very clean film of a Pippard superconductor (In with mean free 
path 11 /im) was investigated using simultaneous magneto-optical imaging and electrical transport 
measurements. A consistent framework for the IS in a flat slab is presented for the first time. 
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The rich physics of the intermediate state (IS) in type- 
I superconductors and unsolved fundamental problems 
continue to keep the IS in a focus of interests of super- 
conductivity P, 0] ■ Two of those problems are the mag- 
netic flux density B in the normal (N) domains and the 
critical field Hd for the IS-N transition. On the other 
hand, the IS provides access to one of the most funda- 
mental parameters, i.e. the Pippard coherence length 
(size of Cooper pairs) , however an experimentally verified 
recipe to extract from the IS parameters is still miss- 
ing. Similar problems, i.e. B in vortex cores, or size of 
the cores, and extraction of microscopic parameters from 
parameters of the mixed state (MS), are among central 
problems of unconventional superconductivity Q. Due 
to an ultimate relationship between the IS and the MS 
one can expect that solution of the problems for the 
former will help to solve them for the latter. 

Ever since Landau introduced a laminar model (LLM) 
for a slab in a transverse field Q , many theoretical mod- 
els of the IS have been proposed, however no one of those 
is adequate. Here we report on an experimental study of 
the IS resulting in the development of a new model for 
the slab in an arbitrary oriented magnetic field. The 
model, based on Tinkham's interpretation of the LLM 
0, is surprisingly simple. Nevertheless it correctly ac- 
counts for the main properties of the IS including ^q- 

When a type-I superconductor with demagnetizing fac- 
tor T] 6] is subjected to a weak magnetic field H, the 
sample is in the Meissner state until H reaches Hi = 
(1 —r])Hc, He being the thermodynamic critical field. At 
Hi the sample undergoes a transition to the IS, where it 
brakes up into N and superconducting (S) domains with 
flux densities B and zero, respectively. Under increasing 
H the normal fraction pn = Vn/V {Vn and V are the 
volumes of N domains and of the sample, respectively) 
increases until the entire sample becomes normal at Hd- 
Usually one assumes that Hd is equal to or slightly 
less than [8j] He- Below it is shown that this is true only 
for very thick samples. 

The shape of the IS domains depends on a big variety 
of factors (sample geometry, material, purity, etc) and is 
unpredictable due to small energy difference for different 
domain shapes However, the shape is simple in two 



cases: a long cylinder in the transverse field and a flat 
slab in a strongly slanted fleld. In the former ?7=0.5 and 
domains are discs perpendicular to the sample axis. In 
the latter ri=l (Hi=0) and domains are linear laminae 

[1) . IS in the flat slab can be investigated using magneto- 
optics (MO) [9|; this is an experiment we performed. 

Primarily we focus on the following questions. (1) 
How does the flux density B and the critical field Hd 
depend on parameters of material, and how do these 
quantities evolve with magnitude and orientation of H? 

(2) How can the domain-wall energy parameter S (and 
therefore also the Ginzburg-Landau (GL) ^(T) and the 
Pippard coherence lengths) be inferred from the lam- 
inar structure? Relationship between S, £_{T) and de- 
pends on material and its purity. For Pippard super- 
conductors (materials with the GL parameter k 1) 
in the pure limit (the elastic mean free path ^ ^ £q) 
(5(r) = 1.89C(r) = L40Co/(l - 0° •^ where t = T/T^ Q. 

On the first glance the answers are available in classic 
texts 0, [^-[1] . However it turns out that in In, directly 
determined from combined Low-Energy Muon Spin Ro- 
tation (LE-/iSR) and Polarized Neutron Reflectometry 
(PNR) experiments (380 nm) [l^l, differs from that fol- 
lowing from d obtained from the IS (240 nm) [llj]. This 
signals for inadequacy of theoretical models used for in- 
terpretation either of the LE-/iSR/PNR data or of the 
IS data. This discrepancy is even more intriguing if one 
takes into account that the IS studies [HI, [12 "played an 
important role in establishing the GL theory in the first 
place" 0. To resolve this contradiction was an original 
motivation for this work. 

The IS structure in the fiat slab was first treated by 
Landau in 1937 [H|- In this work Landau established the 
concept of a surface tension of a S/N interface 7 {5 was 
later defined as Sttj/H^ [l^) and proposed the LLM. 
Assuming that B ^ He, Landau calculated the shape 
of rounded corners of a cross section of the S-laminae. 
The rounded corners yield an excess energy of the system 
competing with the S/N interface energy. Minimizing a 
sum of these two energy contributions. Landau obtained 
the period of the laminar structure D ~ \/'^dJJiXfin) , 
where d is the sample thickness and fL{Pn) is a spacing 
function with pn = h = H/Hc [21 • 
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Soon thereafter Landau, advised by Peierls, admitted 
that LLM is unstable because roundness of the corners 
makes B near the surface less than He- Landau then 
introduced a "branching" model flS*!, which was shown 
to be inadequate by Meshkovsky and Shalnikov in 1947. 
Due to that in 1951 Lifshitz and Sharvin returned back to 
LLM and calculated numerically fL{Pn) and B; identical 
results have later been obtained analytically by Fortini 
and Paumier [7]. i? in LLM, being He in the bulk, near 
the surface is 0.66i?c at low H and increases with H up 
to He at the IS-N transition. 50 years later direct bulk 
/iSR measurements of i? in a Sn slab in a transverse field 
revealed that B{H) starts from He and decreases with H 
down to Hci [14]. Such a dependence for B{H) had been 
anticipated by Tinkham 

De Gennes [Sj] noticed that a positive 7 should reduce 
Hci- Assuming a small reduction, De Gennes obtained for 
the transverse configuration Hci± = Hc[l — 0.9{6/d)'^'^]. 

Tinkham [6] modified LLM in a different way. He no- 
ticed that the dominant contribution to the excess en- 
ergy term comes from inhomogeneity of the field outside 
the sample near its surface and therefore contribution 
due to rounded corners can be neglected. Tinkham com- 
puted this energy introducing a "healing length" Lkeai 
over which the field relaxes to its uniform state away 
from the sample: L'f^eai ~ + ^7^: where -D„ and Ds 
are the width of the S- and N-laminae, respectively. For 
the same reason Tinkham assumed a rectangular cross 
section of the laminae and hence a uniform B, which was 
allowed to be < He- Tinkham's expression for D has a 
form similar to that in LLM with fhiPn) substituted by 
p^(l — Pn)^. For Hci Tinkham obtained 

He^l_ = He[{l + ^S/df'^ ~ 2{5ldf-^] (1) 

with no preliminary restrictions for the Hd reduction. 

The structure expected from the LLM has never been 
observed. Images, reported from the 1950s on revealed 
far more complicated flux patterns than that in the LLM 
and that domains are often corrugated laminae [9] . Bal- 
ashova and Sharvin in 1956 and Faber ^15i] proposed that 
the corrugations stabilize the laminar structure. However 
the disordered laminae in the perpendicular field might 
imply that asymmetry in the field configuration is able to 
change the structure in favor of the Landau's prediction, 
as was first predicted by Dzyaloshinskii in 1955 p^. 

Sharvin (12j was the first to observe the ordered lami- 
nar pattern in a strongly tilted field on 50 mm x 2 mm (di- 
ameter X thickness) Sn [T2l| and In fll'] samples. From the 
images, obtained utilizing a Bitter method 9], Sharvin 
measured D at different temperatures and calculated 
5{T) using an extended LLM, which assumes B\\ — H\\ 
{B\\ and H\\ are in-plane components of B and 7J, re- 
spectively) on top of Landau's original assumption that 
B = He- Sharvin's equation for D has a form 

D^ = J^ ^ (2) 




FIG. 1: MO images taken at 2.5 K. Superconducting regions 
are black. in Oe, H± in Oe): a, (0, 1); b, (60, 8); c, (100, 
6); d, (110, 3); e, (115, 1.3). 



The value of S{0) obtained for In using Eq. (|^ flll| is 
that from which we started our story. Faber [15| crit- 
icized Eq. ([2]) arguing that may alter the shape of 
the laminae corners. One can add that if = the 
magnetic flux is not conserved and therefore the energy 
balance in the system has to be reconsidered. 

From what was said follows that validity of the ex- 
tended LLM [l3] is questionable. Therefore the value of 
^0 following from (5(0) obtained in Ref.O is questionable 
as well. On the other side, the 6{T) dependence obtained 
by Sharvin is correct because it agrees with the GL the- 
ory. And on the third side, the question how to extract 

6 from the IS pattern in the tilted field remains to be 
answered. The last two issues along with the questions 
on B and Hd are addressed below. 

The MO imaging was achieved using a set-up equipped 
with a 3D vector magnet. The sample was a 2.5 /xm thick 
In film deposited on a Si02 wafer. The film residual re- 
sistivity ratio was 540, and the elastic mean free path 
£ w 1 1 fim. The rest film characteristics were the same 
as those reported in Ref.fl^. Overall, the film was a Pip- 
pard superconductor (k ~ 0.07) in the pure limit. The 
sample length was 1 mm and the ratio width/thickness 
was 120, meaning that for the transverse and longitudinal 
fields 7y is 1 and 0, respectively. Images were taken simul- 
taneously with measurements of the electrical resistance 
R using a small low frequency (11 Hz) AC current. 

Typical images are presented in Fig. 1. Contrary to Pb 
, no flux tubes were observed in our In sample both for 
increasing and decreasing H±. At H± < 1 Oe fraction- 
ated laminae were seen in some runs. In the transverse 
field the pattern was disordered laminae (Fig. la) inde- 
pendent of the history of the applied field. This suggests 
that the laminar structure is a ground state topology of 
the IS. The same conclusion was drawn for Al [l^. Dif- 
ference in the pattern in Pb from the pattern in In and 
Al can be attributed to a large difference in k (0.6 for Pb 
and 0.01 for Al) and, accordingly, in 7 for these metals: 

7 in Pb is much smaller than that in Al and In Q . 
While Dn (Ds) increased (decreased) and R varied lin- 
early with H±, the period D = Dn + Ds, being depended 
on was constant at 0.3 < H±/Hic± < 0.8. Near 
H± = and Hci± number of the laminae H decreases. 
The IS-N transition at decreasing field was accompanied 
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by supercooling of the N-state. This confirms that the 
IS-N transition is a first order phase transition 3 . 

Fig. 2 presents the sample phase diagram Hc{T) — 
H^\f^{T) measm'ed with a DC magnetometer. It is com- 
pared with data we measured on other samples and with 
data from Ref.llTl The Hc{T) curve was used for in situ 
determination of the sample temperature. The lower 
curve presents Hci±{T) at = 0, determined on dis- 
appearance of the last S-lamina in MO images and from 
R{H±) measurements; the two perfectly coincide. The 
stars represent Hci± calculated from Eq. ([T]) with — 
380 nm and are clearly consisted with the experimental 
data. De Gennes' formula yields values of Hci± consider- 
ably exceeding the experimental data. Hence Tinkhani's 
interpretation of the LLM is used for further analysis. 

The model. Since we adopt Tinkham's approach and 
in agreement with the experimental images, the N- and 
S-domains are assumed to be rectangular parallelepipeds 
extending in the direction. The contribution of the 
negative surface tension at the S / vacuum interface is ne- 
glected since the penetration depth X <^ d. The out-of- 
plane and in-plane demagnetizing factors are ri±=l and 
?7||=0, respectively. The former means that B±pn H± 
(conservation of the flux of the out-of-plane component 
of the magnetic field), whereas the latter means that 



H\\ and therefore the flux of Bii is not con- 



served. Hence an appropriate thermodynamic potential 
is F ^ F - T/(S||iJ||/47r) = F - Vn^Hj/AT:), where F is 

free energy Q. Summing up the sample free energy at 
zero field, energy of the field B in the sample, energy of 
S/N interfaces and the excess energy of the field over the 
healing length, one obtains for / = F /V: 
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where fno is the free energy density of the sample in the 
normal state at zero field. 

Similar to the case of the LLM, competition between 
the last two terms provides the equilibrium D: 
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Substituting in ([3]) one obtains 
Hi , . 
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f = fnO + ^ [Pn - 1 



(4) 



Pnhi + Ah^J -{1 - p,,)U5) 



where h± = H±/Hc and = H\\/Hc. 

Eq. ([5]) has a minimum yielding equilibrium p„ 

pl = h\/{l-Ah^^d-h\). 
At SEN transition hence 



hc^l_ = ,jA(J/d) + l^^^ - 2^/5jd. 



(6) 
(7) 
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single crystal 
toil 0.5 mm 
Eq. (1)witti 
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FIG. 2: The phase diagram of our sample in the parallel and 
perpendicular fields. The star symbols correspond to Hci_L 
calculated from Eq. ([T|) using ^o=380 nm. The full curves are 
parabolic fits of the data obtained with our sample. 



Finally, the reduced fiux density b ~ B / He is 



(8) 



Obtained expressions satisfy limiting cases, i.e. p„ — > 



at H± 0; Hci_ 
pies i^/SJd < 1) 



_ — > at He- In very thick sam- 

B ^ He and Eq. (g]) converts to Eq. @ 
if — p„)^ is replaced by fL{Pn)', this explains correct- 
ness of the temperature dependence S{T) in Refs.HTI, [121. 

For a perpendicular field (iJy — 0) we have that (a) 
according to Eq. ([5]) , B decreases with increasing H from 
He down to Hci, in agreement with the p,SR results p^ : 
and that (b) Eq. ([7]) reduces to Eq. H]), implying that the 
theoretical points in Fig. 2 are the same in our model. 
Hence our model, composed for the regular laminar struc- 
ture, can be used for irregular laminar patterns as well. 

Interestingly, we find that in a tilted field B is deter- 
mined by H± (Eq. ([5])). This implies that when H± is 
fixed and H\\ is varied, pn (the fraction Dn/D) in Eq. ^ 
varies in such a way that the total B remains constant. 
In Fig. 3 the experimental data for Hci± at nonzero 
are compared to the Hci±{H\\) dependence given by 
Eq. ([7]). We find that Eq. ([7]) correctly describes the ex- 
perimental data. 

Fig. 4 presents the results for 5 at T=1.7K obtained 
in three ways. The circles represent S calculated from 
Eq. ([2]) and Eq. ([4]) using our experimental data for aver- 
age Z? = 1/H and H± at p„ = 0.5. A dashed line in Fig. 4 
represents S directly calculated using — 380 nm. We 
find that the results following from Eq. (g]) agree with the 
directly calculated value of 6. However it is not the case 
for 5 obtained using Eq. ^ . 

Fig. 5 presents experimental data for the average D 
obtained in two runs at T=1.7K and — 0.85i/c and 
corresponding theoretical expectations following from 
Eqs. and Q. Solid circles are the data obtained at 
increasing H± with the sample cooled at H± = 0. The 
data represented by open circles were taken after H± was 
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eq. (7) with 
^o=380nm 
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FIG. 3: The transverse critical field He 
= HJHc Hei±{0) is Hci± at Hn = 0. 
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FIG. 4: The wall energy parameter 5 inferred from MO images 
taken at r = 1.7K using Eq. (01 and Eq. (g)) plotted on two 
scales. The dashed line represents the directly calculated S — 
1.4Co/(l - T/Tcf-^, Tc = 3.415 K and = 380 nm, with the 
error bar corresponding to the error of flQ^ . 



reduced from above Hci± down to 0.5Hci±- In Eq. ^ D 
is controlled by /l; in our model the spacing function is 
f = {I - pl)h\ with pn and h±_ linked by Eq. ©. We 
find that both and / are qualitatively consisted with 
the experimental data, however " works" better at the 
high end of the IS, whereas / is better at the lower end. 
Therefore a refined model should also include contribu- 
tion of the rounded corners. 

Summarizing, performed magneto-optical study of the 
IS in the pure limit Pippard superconductor resulted in 
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FIG. 5: Experimental data (circles) and theoretical curves for 
the period of the laminar flux structure. 



the development of a new theoretical model of the IS for 
a flat slab in an arbitrary oriented magnetic field. The 
model reproduces experimental data on the critical field 
Hci for transverse and tilted field directions; the directly 
measured dependence B{H±) [Ij]; and the directly mea- 
sured value of the Pippard coherence length The 
model allows one to infer S from IS patterns in a strongly 
tilted field, where laminae are linear, and to make a plau- 
sible estimate of 6 from Hd at the transverse field, where 
laminae are corrugated. Overall, the model is a good first 
order approximation of the IS in a flat slab put forward 
by Landau more than 7 decades ago. 

A weak point of our model, which should be addressed 
in a refined approach, is an oversimplified form of Lheai- 
This is a major reason of discrepancy between experimen- 
tal and theoretical D(h±) (Fig. 5), deviation of Pn{h±) in 
Eq.(6) from a linear dependence following from linearity 
of R{H±), and deviation of (5(/ij|) from a constant value 
at large (Fig. 4). To resolve this issue measurements 
of the magnetic field near the surface outside and inside 
the sample are desirable. We note that this is also rele- 
vant for elucidation of the vortex core structure in type-II 
superconductors [18|. 
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